Multielectron dissociative ionization of molecules by intense laser radiation 



Miroslaw Brewczyk 
Filia Uniwersytetu Warszawskiego, 
ul. Lipoma 41, 15-424 Bialystok, Poland 



> 

o 

(N 



O ■ 



> 
o 



On 
O 

Oh' 



X 



Kazimicrz Rzazcwski 
Centrum Fizyki Teoretycznej PAN and College of Science, 
Al. Lotmkow 32/46, 02-668 Warsaw, Poland 

Charles W. Clark 
Electron and Optical Physics Division, Physics Laboratory 
National Institute of Standards and Technology, Technology Administration, 
U.S. Department of Commerce, Gaithersburg, MD 20899 



We solve the hydrodynamic-ballistic equations of motion 
for a one-dimensional time-dependent Thomas-Fermi model 
of CI2 exposed to an intense subpicosecond laser field, and 
observe simultaneous multielectron ionization and molecular 
dissociation. The fragment kinetic energy defect with respect 
to the simple Coulomb explosion picture is found to originate 
in ejected-electron screening of the escaping fragments; its 
magnitude agrees with that observed in recent experiments. 

PACS Numbers: 31.20.Lr,33.80.Gj,33.80.Rv 



Molecules exposed to intense (/ > 10 13 W cm -2 ) sub- 
picosecond laser pulses undergo multiple ionization ac- 
companied by dissociation, a process denoted by multi- 
electron dissociative ionization (MEDI) jl] || . The sim- 
plest picture of MEDI is that of the "Coulomb explo- 
sion," in which the laser field quickly strips a number of 
electrons from the molecule, and the charged ionic frag- 
ments then dissociate under the influence of the repulsive 
Coulomb potential. In this picture, the total kinetic en- 
ergy of the fragment pair with charges q\e, q±e should 
be equal to T — qiq2e 2 /(4TreoR e ) (SI units), where q\ 
and (J2 are integers, e is the elementary charge, and R e 
is the equilibrium internuclear separation of the neutral 
molecule. 

In experiments, however, it is found that the fragment 
kinetic energies are consistently lower than those pre- 
dicted by the Coulomb explosion model, by an amount 
called the kinetic energy defect, A. This appears to be 
true for all observed fragmentation channels, and does 
not depend much on other experimental parameters such 
as laser wavelength and pulse duration. The fragment 
energies are instead consistent with a Coulomb explosion 
at a value of R some 20-50% greater than R e . However, 
since the characteristic time scale for molecular internu- 
clear motion (a vibrational period) is long compared to 
the time scale on which strong-field multiple ionization 
occurs, it is difficult to see how a simple Coulomb explo- 
sion picture could be applicable. 

This Letter puts forth a simple explanation of these 
MEDI observations, and supports it by numerical calcu- 



lations applied to a model molecule. A Coulomb explo- 
sion is initiated at R e , when the most loosely-bound elec- 
trons are stripped from the molecule, but the ejected elec- 
tron cloud does not expand rapidly. Ionization continues 
until only a tightly-bound electronic core remains, and is 
then shut off by the strong internal molecular field (al- 
though for very intense radiation there may be some sub- 
sequent ionization of the separated fragments, i.e. post- 
dissociation ionization (PDI)). Most electron stripping 
thus occurs near R ss R e , so one sees a brief interval of 
acceleration of the ionic fragments by increasing mutual 
Coulomb repulsion. However, when the ionization shuts 
off, a countervailing decelerating tendency becomes ap- 
parent: as the fragments separate through the ejected 
electron cloud, more electronic charge is encompassed by 
the increasing internuclear separation R(t), and the resul- 
tant screening slows the fragments. This post-explosion 
screening gives a value of A consistent with that observed 
in experiments. Charge-symmetric fragmentation is con- 
firmed Q to be a major channel of dissociation. These 
conclusions differ from those of other models that have 
been proposed recently [^],[|-[6j . 

We have done calculations to illustrate this picture, 
by solving the time-dependent equations of motion for 
the nuclei and electrons of a model CI2 molecule in a 
strong laser field. All electrons are treated explicitly, via 
time-dependent density functional theory. We know of no 
previous treatment of molecules in strong fields that deals 
explicitly with more than one electron. Our approach is 
made as follows. 

In strong-field laser irradiation of diatomic molecules, 
it is believed ^| that the molecular axis is quickly aligned 
with the polarization vector of the radiation field, so that 
the ionization process can be modelled with the molecu- 
lar axis parallel to the electric field. To simplify the nu- 
merical calculations, we make the additional approxima- 
tion of confining all particles in the system to move in one 
dimension, defined by this axis. This approximation has 
been used to model strong-field interaction with one-and 
two-electron atoms [jjj and with the H^" molecule || , and 
it has replicated much of the essential physics of three- 
dimensional systems. We utilize time-dependent den- 
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sity functional theory to deal directly with the electron 
density vs. the many-electron Schrodinger wavefunction. 
We have previously applied time-dependent density func- 
tional theory to treat three-dimensional atoms in strong 
radiation fields M , and obtained results for multiple elec- 
tron ionization that agree well with experiments. 

The electron density of the molecule is described as a 
fluid of mass density p(x,t), with an associated velocity 
field v(x, t), which obey the equations of motion: 



dv dv 19 e d 
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The first of the eqs. (|l|) is the usual continuity equa- 
tion. The second is the classical equation of motion for 
an infinitesimal element of fluid subject to forces due to 
the gradients of an electrostatic potential $ and a pres- 
sure P — —(dU/dx) S) where U is the internal energy 
density of the fluid. The effects of quantum mechanics in 
this equation must be expressed by a constitutive relation 
between P and p. Such a relation was derived in a semi- 
classical approximation by Thomas and Fermi (TF) jl0| ; 
its variants have been used to investigate the ground state 
properties of atoms Jll],[l2j and weak-field atomic radia- 
tive response |0|. The TF model treats the electrons 
as a Fermi gas at temperature T = 0, and determines 
the energy of the gas by filling the available phase-space 
volume, subject to the Pauli exclusion principle. This 
gives U of the three-dimensional electron gas as a local 
function of the density, of the form U(r) ~ p(r) 5 / 3 . 

We have sought a corresponding relationship applica- 
ble to the one-dimensional many-electron system. The 
Hohenberg-Kohn theorem (l4j ensures that the energy 
of a stationary one-dimensional system is a functional of 
the electron density, but the semiclassical TF arguments 
may not provide a good approximation to the true en- 
ergy functional. Experimentation with possible forms of 
the functional leads us to propose the following expres- 
sion for the energy density of a one-dimensional electron 
gas, U(x) = An 2 (x), where n{x) is the electron density 
function (n(x) — p(x)/m), and A is a universal constant. 
The value of A is determined by requiring that it give 
reasonable physical properties (binding energies, etc.) of 
the model systems. Such an expression is also obtained 
by applying the semiclassical method in two dimensions 
& 

In treating Coulomb interactions in one dimension, it 
is necessary to eliminate the singularity in the potential 
at x = 0. We follow a standard procedure and 
write the electrostatic potential <J> of our model diatomic 
molecule as 
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Here Z is the atomic number and x\{t),X2{t) are po- 
sitions of the two nuclei. The parameters b, c serve to 
smooth the Coulomb interaction at points of coalescence. 

By solving the static case of eq. |J we obtain the lin- 
ear relationship, n[x) — e<&(x)/(2A), which provides an 
integral equation for the electrostatic potential of the 
time-independent solution Jl6| . Its linearity allows us 
to construct the solution for a many-atom system by lin- 
ear combination of appropriately shifted single-atom so- 
lutions. It is easy to show that the parameters A, b, and c 
are not independent pq |: A = e 2 ln(c/6). Thus we always 
have c > b. Furthermore, our numerical calculations of 
values of R e and D for various values of Z indicate that 
the value of A must be essentially independent of Z, and 
that A = 0.3 a.u. [in the usual atomic units (a.u.) in 
which the numerical values of the electron mass m, the 
elementary charge e, and reduced Planck's constant h are 
equal to unity] . If A lies outside a narrow band of values 
around 0.3 a.u., one obtains unreasonable values of R e 
and D. 

We have applied this model to describe the CI2 (Z = 
17) molecule in a strong radiation field. By solving the 
static equation as a function of internuclear distance, we 
find the Born-Oppenheimer potential curves for CI2 and 
CI2" displayed in Fig. 1. These results were obtained 
with the parameter set (A = 0.3 a.u., b = 1.22 a.u., 
c = 1.65 a.u.) which yields R e = 3.8 a.u. and D = 0.11 
a.u., in agreement with the experimental values fl7j] of 
i? e = 3.76 a.u. and D = 0.091 a.u. respectively. 

We treat the time-dependent system by adding the 
electric dipole interaction with the radiation field to (g) , 
and solving the hydrodynamic equations for the elec- 
tronic density simultaneously with the classical equations 
of motion for the positions of the nuclei. We treat the 
equations (Q) as an initial value problem, with the density 
p(x, t — 0) given by the static solution, and the nuclei ini- 
tially at rest and separated by R e . The first check of the 
time-dependent method is to verify that it reproduces 
molecular vibrations when no external field is present. 
That it does is shown in the inset of Fig. |l|, which dis- 
plays the free oscillation of the internuclear distance in 
the absence of a radiation field. The period of oscilla- 
tion of 34 fs is shorter than the experimental value Jl]]] 
of 59 fs, so our model overestimates the stiffness of the 
molecular potential at R e . 

We now discuss the effects of an intense subpicosecond 
laser pulse on this system. The laser frequency of u> = 
0.0746 a.u. (wavelength A = 610 nm), is that used in 
the experiments of refs. [^|J|]. The laser pulse is turned 
on with a sin 2 (7rf/(2r)) ramp; we treat short-pulse and 
long-pulse cases in which the ramp time r is 10 or 70 
optical periods T (20 or 140 fs) respectively. 

Fig. ^ shows short-pulse results for a peak field 
strength F p — 0.6 a.u. (corresponding to a laser intensity 
I = 1.3 x 10 16 W cm" 2 ). The inset to Fig. § (a) shows 
the Coulomb explosion: the nuclei remain essentially at 
rest until a time just before F p is reached; then there is 
rapid acceleration (increase in cycle-averaged kinetic en- 
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ergy) during an interval of about 10 T; and then a gradual 
deceleration. The dynamics of this system can be read- 
ily understood in terms of the time dependence of the 
charge distribution, which is summarized in Fig. |^(b). 
The lowest, nearly straight line, labelled "grid" displays 
the total electric charge on our grid, and thus represents 
the number of electrons that have been absorbed at the 
boundaries; the "atomic ion" value consists of the charge 
within 3 a.u. of either of the ionic fragments, which is 
seen to be a good long-time indicator of the charge lo- 
calized on that fragment; and "molecular ion" indicates 
the sum of the two atomic ion charges, plus all charge 
in between. We can see that by the time the molecu- 
lar and atomic ions have settled into steady behavior, at 
t w 20T, no charge has been absorbed at the walls of our 
box; the gradual departure of charge from the system at 
later times does not significantly affect them. Note that 
the final charge state of each atom is +4, whereas that of 
the "molecule" is +5. Thus there are three units of neg- 
ative charge (electrons) left in the space between the two 
fragments, which screen their mutual Coulomb repulsion. 

The molecular charge is seen to show a sudden de- 
crease at t » 10T, which corresponds to an increase in 
the enclosed number of electrons. This results from two 
effects: the shutoff of molecular ionization; and the in- 
creasing distance between the fragments, so that more of 
the ejected electron cloud lies between them. It is not 
due to a rescattering mechanism, which would occur on 
the time scale of 1 field period. The mechanism of ioniza- 
tion shutoff is shown in Fig. [|, which displays snapshots 
of the self-consistent molecular potential e<J>(a;,i) early 
in the pulse. The two atomic potential wells are super- 
posed on the linear potential of the electric field, which 
oscillates in time. Electronic charge on the sides of the 
atomic wells or in the interestitial region will be stripped 
off early in the pulse, since there is no barrier to its es- 
cape. However, by t = 8.75T the loosely-bound charge 
is gone, and the remaining charge is stuck in the atomic 
wells (tunneling does not occur in our model; it would be 
slow in real systems under these conditions). The shutoff 
of the ionization thus stops the molecular charge from in- 
creasing; the actual decrease of charge that follows is due 
to the enclosure of previously-ejected electronic charge 
between the separating ion fragments. Comparison of 
frames (a) and (b) of Fig. || shows that roughly two 
units of electronic charge are added to the molecule in 
this way within about 20 fs (10 T) after ionization shut- 
off, at which time the deceleration begins. This results 
in A w 20 eV, in agreement with the experiments of ref. 

For completeness we make some observations on long- 
pulse case. If the ramp of the pulse is increased, we see 
essentially the same deceleration mechanism as in the 
short-pulse case, although PDI occurs at high intensities. 
Fig. ^ shows molecular charges for two intensities, at the 
same frequency as in Fig. |2j, but with a ramp time of 
140 fs. At both intensities we see a knee in the curve 
similar to that in Fig. |[ and it occurs at about the same 



value of instantaneous field strength. However, at the 
higher intensity, the knee represents only a temporary 
disruption of ionization, and PDI occurs as the field rises 
to its peak strength. 

In conclusion, we have presented the first dynamical 
model of MEDI. The results are in qualitative agreement 
with experimental data. Reproducing the main features 
of the experiments, the model does not support the idea 
of stabilization Q , nor does it attribue much importance 
to the enhanced ionization at some critical distance ||. 
We predict some post-dissociation ionization for longer, 
picosecond pulses with peak intensities above 10 16 W 
cm~ 2 . This could be tested experimentally. 
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FIG. 1. Potential energy curves for Cb and ClJ 2 , as com- 
puted with the parameters A = 0.3, b = 1.22, c = 1.65 a.u. 
Inset: Internuclear distance as a function of time for free os- 
cillations of CI2 for small internuclear excursions away from 

Re- 

FIG. 2. Evolution of the molecular features for F p — 0.6 
a.u. and uj = 0.0746 a.u. with the field ramped to maximum 
intensity in 10 optical periods T, as marked by the arrow, 
(a) Kinetic energy of the atomic ion fragments vs. time; in- 
set shows internuclear distance vs. time. Note that Coulomb 
explosion begins near the peak intensity; deceleration of the 
fragements sets in around t — 20T. (b) Distibution of net 
electrical charge (includes both nuclear and electronic com- 
ponents). Grid: net charge of the system on our spatial grid, 
which is contained in a box of length 200 a.u.; increase around 
t = 2QT comes from ejected electrons leaving the box, as sim- 
ulated by absorbing boundary conditions. Atomic ions and 
molecular ions: see text. 

FIG. 3. Self-consistent potential curves for the case of Fig. 
I at t = 4.25, 7.25, and 8.75 T, as indicated. 

FIG. 4. Molecular charge vs. time uj — 0.0746 a.u. in the 
case of a long pulse ramp, at two values of peak intensity as 
labelled: F p is attained at 140 fs as marked by arrow. The 
effect of ionization shutoff and ejected-electron enclosure are 
similar to those seen in the short-pulse case of Fig. § The 
knee occurs at roughly the same instantaneous field strength 
in the two cases, corresponding to the shutoff condition dis- 
played in the bottom frame of Fig. [| PDI is seen at the 
higher intensity 
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